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We study transport through a quantum dot coupled to normal and superconducting leads 
using the numerical renormalization group method. We show that the low-energy properties of 
the system are described by the local Fermi liquid theory despite of the superconducting corre- 
lations penetrated into the dot due to a proximity effect. We calculate the linear conductance 
due to the Andreev reflection in the presence of the Coulomb interaction. It is demonstrated 
that the maximum structure appearing in the conductance clearly characterizes a crossover 
between two distinct spin-singlet ground states, i.e. the superconducting singlet state and the 
Kondo singlet state. It is further elucidated that the gate-voltage dependence of the conduc- 
tance shows different behavior in the superconducting singlet region from that in the Kondo 
singlet region. 
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1. Introduction 

Recent progress in nano-technology has attracted 
much interest in studying quantum transport in meso- 
scopic systems. Among others, a quantum dot (QD) 1 has 
played an important role to reveal correlation effects in 
the nanoscale systems. In particular, the observation 2,3 
of the Kondo effect in QD systems 4-8 opened a sys- 
tematic way to investigate strongly correlated electrons, 
which has encouraged further theoretical and experimen- 
tal studies in this field. Besides this substantial progress, 
transport properties of a mesoscopic system with hybrid 
normal(metal)-superconductor junctions have been also 
investigated extensively. In this system, the Andreev re- 
flection 9 plays a key role in physics, in which an incom- 
ing electron from normal side can be reflected as a hole, 
therefore transferring a Cooper pair into superconductor. 

The above interesting topics in nanoscale systems nat- 
urally stimulated the research on the Andreev reflec- 
tion for a QD coupled to normal and superconducting 
leads (N-QD-S). 10-15 In this system, the Andreev reflec- 
tion (the proximity effect) at the QD-S interface induces 
the superconducting correlation in the QD, which has 
a tendency to form the spin-singlet state. On the other 
hand, for large Coulomb interactions, the Kondo effect 
is enhanced and therefore the Kondo spin-singlet state 
is stabilized between the spin moment in the QD and 
the conduction electrons in the leads. Thus, the compe- 
tition between these two distinct spin-singlet states oc- 
curs in the N-QD-S system. In order to clarify how this 
competition affects the transport in this system, theoret- 
ical analyses 16-27 as well as experimental investigations 28 
have been done intensively, e.g. the linear 16, 17:19,21:25,28 
or nonlinear 16, 18,22,23,25,27,28 conductance, the excess 
Kondo resonance coming from the novel co-tunneling 
process (Andreev- normal co-tunneling), 20 the Andreev 
reflection through a QD embedded in an Aharonov- 
Bohm ring 24 and the adiabatic electron pumping 26 etc. 



In particular, some theoretical studies on the linear con- 
ductance have clarified that the Coulomb interaction 
suppresses the Andreev reflection at the QD-S inter- 
face, which leads to the decrease of the linear conduc- 
tance. 16, 17, 19, 25 However, these studies have been done 
on the assumption that the Coulomb interaction in the 
QD is sufficiently large. On the other hand, Cuevas et al. 
analyzed the conductance over the range from the non- 
interacting case to the strong-correlation limit, using a 
modified second-order perturbation theory. 21 They found 
that the increase of the coupling between the QD and 
the superconducting lead makes it possible to restore the 
conductance possibly up to the maximum value 4e /h, 
although physical implications of the conductance maxi- 
mum in the presence of the Coulomb interaction were not 
discussed in detail. In an N-QD-S system, moreover, the 
total number of electrons is not conserved, because of the 
superconducting correlation. In such a system, it is not 
obvious whether the low-energy properties are described 
by the local Fermi liquid theory. 

In this paper, we theoretically investigate the trans- 
port in an N-QD-S system with the use of the numeri- 
cal renormalization group (NRG) method, 29,30 which has 
been applied successfully to a Josephson current through 
a QD. 31-33 Applying the Bogoliubov transformation, we 
first show that the low-energy properties of an N-QD-S 
system are described by the local Fermi liquid theory. Us- 
ing the NRG method, we calculate the conductance due 
to the Andreev reflection with high accuracy and thus 
confirm Cuevas et a/.'s results. To understand the be- 
havior of the conductance, the renormalized parameters, 
which characterize the Andreev bound states around the 
Fermi energy, are calculated. From the analysis of the 
ground state properties, we demonstrate that the con- 
ductance maximum clearly characterizes a crossover be- 
tween the superconducting singlet state and the Kondo 
singlet state. 
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This paper is organized as follows. In the next sec- 
tion, we introduce the model and describe the low-energy 
properties in terms of the local Fermi liquid theory. Then 
in §3, we calculate the conductance and the renormalized 
parameters. We discuss how the interplay between the 
Kondo effect and the superconducting correlation is re- 
flected in the transport, by systematically changing the 
Coulomb interaction, the tunneling amplitude and the 
gate voltage. A brief summary is given in the last sec- 
tion. 

2. Model and Formulation 

2.1 Model 

The Hamiltonian of a QD coupled to normal (N) and 
superconducting (S) leads is given by 

H = H d + H d + Hs + -Hjv + Hts + Htn, (1) 

where H d +H d and H$(n) represent the QD part and the 
superconducting (normal) lead part, respectively. Hts 
and Htn are the mixing terms between the QD and the 
leads. The explicit form of each part reads 

H° d = (e d + |) (n d 1), H u d = - l) 2 , 

q,a q 

h n = J2 £kC L c k«> h tn = -m=^° d ° + H - c -)' 



k.(T 



Hts = 



V s 



(U+H.c). 



(2) 



The operator d\ creates an electron with energy e d and 
spin a at the QD, and n d = X^o- d-t^a- Following ref. 29, 
we here write down the QD part H d + H d in such a way 
that the energy of the one-electron occupied state (n d — 
1) is zero. Note that H d +H d = e d n d + Un d ^n d i+const., 
because (n d — l) 2 = 2n d ^n d i —n d + l. In this representa- 
tion, H d gives the energy shift due to the deviation from 
the electron-hole symmetric case (e^ + U/2 = 0). In the 
Hamiltonian for leads, sJ cr (cj i(T ) is the creation operator 
of an electron with the energy £ q {£k) in the supercon- 
ducting (normal) lead. In Hts (Htn), Vs(Vn) is the tun- 
neling amplitude between the QD and the superconduct- 
ing (normal) lead, and Ms(Mn) is the number of lattice 
sites in the superconducting (normal) lead. We assume 
that the superconducting lead is well described by the 
BCS theory with a superconducting gap A = |A|e'^ s , 
where <j>s is the phase of the superconducting gap. 

In what follows, we consider the limiting case of | A| — > 
oo. The essential physics of the Andreev reflection, which 
occurs inside the superconducting gap, is still captured 
in this limit. In this situation, the Hamiltonian (1) can 
be reduced exactly to an effective single- channel Hamil- 
tonian (see Appendix A) 32 - 34 

(3) 



H 



off 



where 



m + m 



Hl C = A, 



sc 



H. 



:d\d\ 



Hn + H- 



TN - 



H.c. 



(4) 



H d c denotes the effective onsite superconducting gap at 
the QD, 



A d = Tse 1 ' 



(5) 



Notice here that the resonance strength between the 
QD and the superconducting lead is given by Ts(e) — 
^^2 q Vg5(e — e q )/Ms, which is reduced to an energy- 
independent constant Ys in the wide band limit. For sim- 
plicity, we set <j>s = in what follows, so that A d (= Tg) 
becomes real. The reduction in the number of the chan- 
nels gives us a practical advantage in the NRG calcu- 
lations, because this method works with high accuracy 
for single-channel systems while the accuracy becomes 
rather worse for multi-channel systems. 

2.2 Effective Anderson Hamiltonian with normal leads 
In the NRG method, the normal lead part is trans- 
formed into a linear chain after carrying out a standard 
procedure of the logarithmic discretization. 29 Then, a se- 
quence of the Hamiltonians is obtained as 



H 



eff 
NRG 



= A^- 1 )/ 2 ( H° d + H u d + Hl c +H n +Htn), 



(6) 



where 



Hn + Htn 

N-l 



tn A / ^ fn+l cr fna + fna fn+1 a 



n=— 1 a 



(7) 

In eq. (7), = d a , and f na for n > is an operator for 
the conduction electron in the normal lead. The hopping 
factor t n is defined by t_i = tiA -1 ' 2 for n = — 1, where 



I2T N DA A 
v = \l , A A 



1 / 1 + 1/A 

2 V 1 - 1/A 

and for the conduction band (n > 0), 



log A, 



tn=D 



1 + 1/A 



1 - 1/A 



n+l 



2 y/l - 1/A 2 «+Vi - 1/A 2 ™+ 3 



(8) 



(9) 



Here, D is the half- width of the conduction band, and r^r 
(= TT^2 k V N S(e— £fe)| £ = £F /M w -) is the resonance strength 
between the QD and the normal lead. The factor A A is in- 
troduced to compare the discretized model with the orig- 
inal Hamiltonian (3) precisely, and it behaves as A A — > 1 
in the continuum limit A — ► 1. 29 > 35 

In principle, we can carry out the NRG calculation for 
the Hamiltonian (6). However, we note that H d c in eq. 
(4), which represents the effective onsite superconduct- 
ing gap at the QD, mixes states with different particle 
numbers. This means that the eigenstates of the Hamil- 
tonian (6) can not be classified in terms of the total num- 
ber of electrons. To avoid this inconvenience, we perform 
the Bogoliubov transformation, 36 which is summarized 
in Appendix B. In the present case the superconducting 
gap is absent ihHn + Htn, so that the Hamiltonian (6) 
can be mapped onto the Anderson model without the 
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onsite superconducting gap, 



N-l 



n— — 1 a 



where 



^ = y(^ + |) . (ii) 

The important point is that the total number of Bogoli- 
ubov quasiparticles J2n=-i So- Inalna conserves. More- 
over, the Hamiltonian (10) is identical to the ordinary 
Anderson model (without superconductivity). Therefore, 
the low-energy properties can be described by the lo- 
cal Fermi liquid theory, even though the original model 
of eq. (6) has the onsite superconducting gap. Equa- 
tion (10) has been obtained on the assumption that 
|A| — > oo. Also for finite |A|, the coupling to the nor- 
mal lead via could make the low-lying energy states 
at \e\ < min(|A|,T/f) be described by the Fermi liq- 
uid, where Tk is the Kondo temperature in the case of 

r s = o. 

There is another point to be mentioned here. By com- 
paring the first term of the Hamiltonian (10) with H® 
in eq. (2), we notice that the parameter (sd + U/2) is 
replaced by Ed- This means that the Hamiltonian (10) 
corresponds to the Anderson model with the energy level 
of the impurity site £~d = Ed — U/2, while the term 
(U/2)(rid — l) 2 due to the Coulomb interaction remains 
unchanged. We discuss the case of Ts ^ (A<j ^ 0) in 
this paper, so that we treat the reduced Hamiltonian (10) 
with Ed > in eq. (11) (so-called asymmetric Anderson 
model). 

2.3 Local Fermi-liquid description 

We now introduce the Green function to formulate the 
density of states (DOS) of the QD. Following eq. (B-ll) 
in Appendix B, we can write down the retarded Green 
function of the QD, after the Bogoliubov transformation, 
as follows, 

G5 TidT ( e ) = «^_ lT ,7-iT( e )+^_ 1 i,7-u( e )- ( 12 ) 
Note that the coherent factors Ud,Vd are real because of 
</>s = 0. Using eq. (12), the DOS of the QD is given by 

Pd{e) = --luvG r d]d] {e) 



«3lmG^_ lti7 _ lT (e)+«2lmG^_ ai7 _ a (e)} 

(13) 

The point we wish to stress is that G^_^ ^ (e) and 
G r ^ 7 _ 1 j( e ) are derived from the generalized Ander- 
son model (10), in which the superconductivity does not 
show up explicitly. 

By using the formula (13), we can describe the An- 
dreev bound states in the QD, which arc induced by the 



Andreev reflection at the QD-S interface. We here focus 
on these states around the Fermi energy (e ~ 0), where 
the self-energy due to the Coulomb interaction E(e) is 
approximately given by E(e) ~ S(0) + e 9S(e)/9e| e=0 . 
Then, the retarded Green function for electrons after the 
Bogoliubov transformation reads 



G 7 _ 



T,7- 



s - E d + iT N - S(e) 



e - E d + iT N 



where 



Ed = z(Ed 



E(0)), T N 
o>S(e) 



zTn, 



z= 1- 



<9e 



e=0 



As the retarded Green function for a hole, G 
is given by G;_ l L ^ a (e) = - (g 



7-ii,7-il v I \ 7-iT,7 

end up with the DOS around the Fermi energy 



7-1J..7- 
lT (- £ ))' 



Pdie) 




+ 



d'-N 



(£-Ed) 2 + r% (£ + E d ) 2 + r% 



(14) 

(15) 
(16) 

, we 
(17) 



We next consider the transport properties in the small 
bias regime, which are governed by the Andreev reflec- 
tion induced inside the superconducting gap. Accord- 
ing to Appendix C, the linear conductance Gv=o — 
dl /dV\v=o at zero temperature is given by 



4e 2 



{Ed/T 



E 



d {i + {E d /v N yy 



(18) 



We see that the conductance is^determined by the ratio of 
the renormalized parameters Ed/T^, which are obtained 
from the eigenvalues of the Hamiltonian (10) at the fixed 
point. 37 

3. Numerical Results 

In this section, we discuss transport properties for the 
N-QD-S system at zero temperature. As mentioned in 
the previous section, we assume that the superconduct- 
ing gap |A| is sufficiently large (|A| — ► 00). In this case 
the excited states in the continuum outside the super- 
conducting gap can be neglected. It describes a situation 
where the superconducting gap is much larger than the 
characteristic energies of the Andreev reflection. 

3. 1 Influence of the Coulomb interaction U 

Let us first discuss how the Coulomb interaction U 
affects the transport properties at zero temperature. To 
this end, we explore the detailed properties of the An- 
dreev bound states, which can be obtained from the 
renormalized parameters computed by means of the 
NRG method with high accuracy. These renormalized 
parameters also determine the conductance, so that we 
can clarify how the conductance is controlled by the An- 
dreev bound states formed around the Fermi energy. 

We observe how the Andreev bound states change their 
characters with the increase of the Coulomb interaction 
U. As discussed in §2, the local DOS of the QD around 
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the Fermi energy is given by eq. (17). In particular, in 
the electron-hole symmetric case (ed + U/2 = 0), Eg = 
r s (= A d ) and u\ = v\ = 1/2 follow from eqs. (5), (11) 
and (B-6). Then, the DOS of the QD, Pdi?), is rewritten 

as 

« (t)=! s{( t -fri' + f* + (£+ f"; + fJ, (i9 » 

where 

f s = z(T s + E(0)). (20) 

It is seen from eq. (19) that Pd(s) in the low-energy region 
is determined by the renormalized parameters Ts and 
Tjv. In the noninteracting case (U = 0), they are reduced 
to the bare ones — IV and Ts — Ts- In this case, the 
Andreev bound states are formed at e = ±r<?, which are 
broadened (finite width IV) by the coupling between the 
QD and the normal lead, as schematically shown in Fig. 
1(a). We shall refer to these Andreev bound states with 



(a) U = (b) U * 




Fig. 1. Andreev resonances around the Fermi energy (Ep = 0) 
in the symmetric case (e^ + U/2 = 0). 



and IV. Figure 2(a) shows the renormalized parameters 
Tjv and Ts as a function of the Coulomb interaction U 
for Ts — r^r, where the parameters T N ( S ) are normalized 
by the bare resonance strength, IV- As shown in Fig. 




5 10 15 20 



finite widths as the Andreev resonances in the following. 
Here we would like to comment on Tj, which gives the 
position of the Andreev resonances shown in Fig. 1(a). 
The Andreev reflection at the QD-S interface gives rise to 
the superconducting correlation in the QD. Taking into 
account eqs. (4) and (5), in the case of |A| — ► oo, we see 
the amplitude of the superconducting correlation in the 
QD is given by the resonance strength Ts(= Ad). Then, 
Ts becomes a parameter indicating the strength of the 
Andreev reflection at the QD-S interface, while it origi- 
nally represents the position of the Andreev resonances in 
the QD. When the Coulomb interaction U is introduced, 
the effective tunneling of electrons between the QD and 
the leads is modified, so that the Andreev resonances 
are renormalized. Namely, the position and the width 
of these resonances, which correspond to the renormal- 
ized parameters Ts and IV respectively, become smaller 
with the increase of U (see Fig. 1(b)). Note that when 
the effect of the Coulomb interaction U is included, the 
strength of the Andreev reflection at the QD-S interface 
is given by Ts instead of Ts- Since we are concerned 
with the DOS around the Fermi energy, the Andreev 
resonances away from the Fermi energy are not shown in 
Fig. 1(b). The overall structure including the high energy 
region can be found in the literature. 16, 19-21,25 

To observe the formation of the Andreev resonances in 
our model, we calculate the renormalized parameters Ts 



Fig. 2. (Color online) (a) Renormalized parameters Ts and Tjv, 
(b) expectation value of the pair correlation in the QD and (c) 
linear conductance as a function of the Coulomb interaction U 
for Tg = Tjv in the symmetric case (e^ + U/2 = 0). Pjv(s) an d 
U are normalized by the bare resonance strength, T^- The NRG 
calculations have been carried out for A = 3.0 and T^/D = 
1.0 X 10 — 3 . Inset of (c): Conductance in the noninteracting case 
(U = 0) as a function of the ratio T$/Yn, where we set = 0. 

2(a), both of Ts and Tjv decrease monotonically with 
the increase of U. The decrease of Tn signals a crossover 
from the charge-fluctuation regime to the Kondo regime, 
as is the case for general N-QD-N systems. On the other 
hand, the decrease of Ts results in the suppression of 
the Andreev reflection at the QD-S interface, implying 
that the superconducting correlation is suppressed in the 
QD. To confirm this implication, we also calculate the 
expectation value of the pair correlation in the QD, which 
is given by 

djdj + ^d^ =--tan- 1 (f s /IV). (21) 

Note that the sign of (djrfj + d^) depends on the def- 
inition of the superconducting gap A in Hs'- namely, 
(djdj +d|(i|) becomes negative because the second term 
in the right hand side of eq. (2) is assumed to be pos- 
itive. As shown in Fig. 2(b), we see that the absolute 
value |(rf|<i| +<2i d+}\ approaches with the increase of U 
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and/or Ts in Fig. 2(a). The Andreev resonances around 
the Fermi energy are thus renormalized with the increase 
of U, as shown in Fig. 1. It is to be noticed that in the 
large U region where the Kondo effect is dominant, these 
two resonances merge around the Fermi energy, forming 
a single sharp Kondo resonance, in accordance with refs. 
16, 19-21 and 25. This is indeed seen from the tendency 
that Ts decreases more rapidly than in Fig. 2(a). 

We next discuss how the conductance changes with the 
increase of U. In the symmetric case (ed + U/2 = 0), the 
conductance of eq. (18) is rewritten as 



G 



V=Q 



4e^ 
~h 



4(r 5 /r 



N) 



{i + (r s /r N rr 



(22) 



The conductance of eq. (22) is a function of the ratio 
^s/^n and has a maximum at Ts/Tn = 1, as pointed 
out in the previous studies. 17,21 Here, we would like 
to mention the relation between the conductance and 
the Andreev resonances shown in Fig. 1. In the case of 
Ts/r^ = 1, which corresponds to the condition for the 
maximum of the conductance, the distance of the An- 
dreev resonances measured from the Fermi energy (Tg) 
is equal to the width of these resonances (Fat)- We also 
consider the case of r s /r w ^ 1. When rs/r^r > 1, the 
distance Tg becomes larger than the width Tjv- Then, the 
DOS of the QD around the Fermi energy gets smaller, 
which results in the^decrease of the conductance. On the 
other hand, when Ts/Tn < 1, the Andreev resonances 
get closer to the Fermi energy, which leads to the en- 
hancement of the DOS around the Fermi energy. At the 
same time, however, this means that the transport due 
to the Andreev reflection is suppressed, so that the con- 
ductance decreases as well as the case of Ts/Fat > 1. 

Figure 2(c) shows the conductance computed as a 
function of the Coulomb interaction U for Ts =Tjy. For 
reference, we also plot the conductance in the nonintcr- 
acting case (U = 0) as a function of the ratio Fs/Fat (in- 
set of Fig. 2(c)), which is given by replacing Ts/Tn in eq. 
(22) with rs/rjy. From the inset, we see that the conduc- 
tance for C/ = has the maximum at Ts /Tn = 1 , where 
the couplings to normal and superconducting leads have 
the same amplitude. In this case, the maximum value of 
the conductance reaches the unitary limit 4e 2 /h because 
of the electron- hole symmetry. With the increase of U, 
the conductance decreases monotonically from 4e 2 //i, as 
shown in Fig. 2(c). This behavior is in agreement with 
Cuevas et aZ.'s results, where they treated a system with 
a finite gap A of the superconducting lead. 21 As stated in 
ref. 21, this is due to the reduction of Ts/Tn, which cor- 
responds to the case of 1^/1^ < 1 stated above. Similar 
behavior of monotonic reduction is observed for smaller 

t s (t s /t n < i). 

Note, however, that somewhat different behavior ap- 
pears in the region of Fs/Fat > 1. Figure 3 shows sim- 
ilar plots of the renormalized parameters and the con- 
ductance for rs/rAr = 5. Let us focus on the small 
U region in Fig. 3(a) and (b). In this region, Ts de- 
creases monotonically, which reduces the absolute value 




U/T 



d,d*)\, as is the case of Ts/Fat = 1. On the other 



N 



Fig. 3. (Color online) (a) Renormalized parameters and Tjv, 
(b) expectation value of the pair correlation in the QD and (c) 
linear conductance as a function of the Coulomb interaction U 
for rs/Tjv = 5. The other parameters are the same as in Fig. 2. 



hand, remains almost unchanged unlike the case of 
Ts/Tn = 1. This result indicates that the position of the 
Andreev resonances approaches the Fermi energy, keep- 
ing its resonance width unchanged. However, when the 
value of Ts gets close to that^of Tn,Tn also begins to 
decrease. For larger U, both Ts and Tn approach 0, as 
is the case of rs/rAr = 1. Comparing Ts and Tn in Fig. 
3(a) with the conductance shown in Fig. 3(c), we see that 
the conductance has the maximum when the condition 
Ts — Tn is satisfied. 

Summarizing the above results, we can say that the 
conductance shows the characteristic [/-dependence ac- 
companied by the maximum structure, around which the 
strength of two effective resonances is exchanged: we have 
Ts > Tn in the small U region (U/T N < 10), while 
Ts < Tn in the large U region (U/T N > 10). In the fol- 
lowing, we demonstrate that the maximum of the con- 
ductance characterizes a crossover between the two dis- 
tinct singlet regions where either the superconducting 
correlation or the Kondo correlation is dominant. 

3. 2 Implications of the conductance maximum 

To discuss the condition that the conductance takes 
the maximum value in more detail, we would like to ob- 
serve how the conductance changes as a function of the 
resonance strength Ts that gives a measure of the su- 
perconducting correlation in the QD, as has been done 
in ref. 21. We find it more instructive to summarize our 
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results in two-dimensional plots in the plane of Ts and 
U, which indeed allows us to elucidate how the electron 
correlations affect the conductance maximum. 

Figure 4 (a) is the color-scale representation of the 
conductance as a function of U /Tm and Ts/Tm- We also 
show the enlarged picture in small U and T s region in 
Fig. 4(b). As shown in Fig. 4(b), the conductance takes 






2 4 6 8 10 



Fig. 4. (Color online) (a) Color-scale representation of the con- 
ductance G/Go, where Go = 4e 2 //i, as a function of U/Tfj and 
Tg/Tff. The dashed line, Fg = U/2, represents the boundary of 
the singlet-doublet transition for Fpf = 0. (b) Enlarged picture 
in the region with small t//Tjv and Fs/^n- 

the maximum at Ts/Tn = 1 when U/Tn = 0, in ac- 
cordance with the inset of Fig. 2. As U/Tm increases, 
the value of Ts/Tn giving the maximum of the con- 
ductance increases linearly along the line of T$ = U/2 
(see Fig. 4). As mentioned above, the maximum of the 
conductance reaches the unitary limit (Go = Ae 2 /h) 
because of the electron-hole symmetry. For instance, 
around U/Tpj — 10 and Ts/Tjy = 5, the conductance 
reaches the unitary limit, as already shown in Fig. 3(c). 

Here we consider the physical implication of Ts = U/2, 
which coincides with the condition that the conductance 
reaches the unitary limit for large values of Ts/Tn and 
U/Tn. To this end, let us examine the limit of — > 0. 
When rjv = 0, the QD is disconnected from the normal 
lead, and only connected to the superconducting lead 
(QD-S). Recall here that the QD-S system is equiva- 
lent to a magnetic impurity model embedded in a su- 
perconductor, which has been studied for dilute mag- 
netic alloys. 38-41 As discussed in refs. 38-41, the ground 
state of the QD-S system is a nonmagnetic singlet or a 



magnetic doublet. Specifically in the limit of |A| — * oo, 
the ground state only depends on the ratio of Ts/U. 
Namely, the ground state is a superconducting spin- 
singlet state for Ts/U > 0.5 or a spin-doublet state for 
Ts/U < 0.5. A transition between these ground states 
occurs at Ts/U — 0.5. 

We now observe what happens for finite Tpf. When 
the coupling between the QD and the normal lead is in- 
troduced, conduction electrons in the normal lead screen 
the free spin moment to form the Kondo singlet state in 
the region of Ts/U < 0.5. Therefore, the ground state 
of the N-QD-S system is always singlet. There are, how- 
ever, two distinct singlets, i.e. one with superconducting- 
singlet character for Ts/U > 0.5 and the other with 
Kondo-singlet character for Ts/U < 0.5. From the above 
discussion, we see that the maximum of the conductance 
in Fig. 4 for large values of rs/r^v and U/Tjq clearly 
characterizes the crossover between these two different 
spin-singlet states. We have drawn this conclusion on the 
assumption that the gap is sufficiently large (|A| — > oo). 
We believe that this conclusion also holds for an N-QD-S 
system with a finite gap A of the superconducting lead, 
though the onset of the transition of the ground state for 
Tn = depends also on A. 32 ' 33 

Before closing this subsection, we display the conduc- 
tance in Fig. 5 in a slightly different way to make the 
comparison with experiments easier; it is plotted as a 
function of the ratio Ts/U for several values of T^/U, 
which are both controllable by changing the voltage be- 
tween the QD and the leads experimentally. 2,3 From 




r s /u 

Fig. 5. (Color online) Conductance as a function of the ratio 
Tg/U for several values of Tjy /U, where the Coulomb interaction 
is fixed as U/D = 2.0 X 10" 2 . 

Fig. 5, we see that as Tn/U gets smaller, the peak 
structure becomes sharper and its position approaches 
Ts/U = 0.5. Thus, the peak structure in the conduc- 
tance for Tn/U = 0.05 and 0.005 in Fig. 5 clearly char- 
acterizes the crossover between the superconducting and 
Kondo spin-singlet states, as discussed above. On the 
other hand, as Tn/U becomes larger, the peak struc- 
ture is somewhat broadened and its position shifts to- 
ward Ts/U > 0.5, as shown for Tjq/U = 0.5 in Fig. 5. 
Note here that the increase of T^/U enhances charge 
fluctuations in the QD, making the Kondo correlation 
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weak. Therefore, the crossover between the supercon- 
ducting and Kondo singlet states becomes smeared, and 
accordingly the position of the conductance maximum 
deviates from Ts/U — 0.5. 

3. 3 Away from the symmetric case 

Finally, we discuss the conductance as a function of the 
energy level in the QD, Ed, away from the electron-hole 
symmetric case (ed+U/2 ^ 0). Note that the energy level 
in the QD can be changed by the gate- voltage control. A 
deviation from the electron-hole symmetric case is rep- 
resented by H® in eq. (2), as stated in §2. In this case, 
the position of the Andreev resonances is given by E# 
in eq. (15) instead of T$- Following the way outlined in 
the symmetric case, we discuss the characteristics of the 
conductance in connection with the ground-state nature 
of the system. 




Fig. 6. (Color online) Color-scale representation of the conduc- 
tance G/Gq (Go = 4e 2 /h) as a function of e^/U and Tg/U. In 
this calculation, we set T N /U = 0.05 and U/D = 2.0 X 10~ 2 . 
The dashed line is {[e d /U + 1/2) 2 + {Ts/U) 2 } 1 / 2 = 1/2, which 
gives the boundary of the singlet-doublet transition for Tat = 0. 

In Fig. 6 we show the conductance in the color- 
scale representation as a function of Ed/U and Ts/U. 
We also draw a half circle given by {(sd/U + 1/2) 2 + 
(Ts/U) 2 } 1 / 2 — 1/2, which denotes the boundary of the 
singlet-doublet transition for Tn = 0. The inside of 
the half circle is the spin-doublet region at Tn = 0, 
which is replaced with the Kondo spin-singlet region in 
the presence of any finite Tjy. On the other hand, the 
outside is the superconducting spin-singlet region. From 
Fig. 6, we see that the conductance for Tg/U < 0.5 
has large values along the half circle. This means that 
the conductance as a function of Ed has two peaks. 21 
Note here that a deviation from the electron-hole sym- 
metric case (£d/U 7^ —0.5) makes the Kondo correla- 
tion weaker. Therefore, the two-peak structure of the 
conductance for Ts/U < 0.5 indicates the crossover of 
superconducting-Kondo-superconducting singlet states. 
On the other hand, for Ts/U > 0.5, the ground state is 
the superconducting singlet state for any values of Ed/U. 
As seen in eqs. (11) and (15), a deviation from the sym- 
metric case drives the position of the Andreev resonances 



(Ed) away from the Fermi energy. In the superconduct- 
ing singlet region (Ts/U > 0.5), Ed/Tjy > 1 is satisfied, 
so that the deviation leads to the reduction of the con- 
ductance, as is the case of F^/Tat > 1 in the symmetric 
case. Thus, the conductance as a function of Ed has a sin- 
gle maximum at Ed/U = —0.5, although this peak is not 
a sign of the crossover between two singlet states unlike 
the peaks for Ts/U < 0.5. 

4. Summary 

We have investigated transport properties of an N-QD- 
S system using the NRG method. Especially, we have 
focused on the limiting case of |A| — > oo with partic- 
ular emphasis on the Andreev reflection which arises 
inside the superconducting gap. Adapting the Bogoli- 
ubov transformation to the simplified model, we have 
first demonstrated that our system with superconduc- 
tivity can be mapped on an effective Anderson impurity 
model without superconductivity, which enables us to 
describe the low-energy properties in terms of the local 
Fermi liquid theory. 

To clarify the influence of the Coulomb interaction on 
the transport due to the Andreev reflection, we have cal- 
culated the conductance as a function of the Coulomb 
interaction U. For the ratio Ts/rjv > 1, the conduc- 
tance has the maximum in its {/-dependence while for 
Ts/rw < 1 it decreases monotonically, which is in accor- 
dance with Cuevas et al.'s results. We have also calcu- 
lated the renormalized parameters to discuss the conduc- 
tance in connection with the Andreev resonances around 
the Fermi energy. Through the analysis using the renor- 
malized parameters, we have found that the maximum of 
the conductance gives an indication of the crossover be- 
tween two distinct types of singlet ground states. To ob- 
serve the nature of the crossover in detail, we have stud- 
ied the transport properties by focusing on the changes 
of the Coulomb interaction U and the resonance strength 
Ts- In particular, starting from the special case with 
= 0, i.e. the QD system only coupled to the su- 
perconducting lead, we have shown that the conduc- 
tance maximum clearly characterizes the crossover be- 
tween the Kondo singlet state and the superconducting 
singlet state. 

It has been further elucidated that the gate-voltage 
dependence of the conductance shows different behavior 
depending on the value of Ts', there are two peaks in the 
gate- voltage dependence characterizing the crossover of 
superconducting-Kondo-superconducting singlet regions 
for Ts/U < 0.5, whereas only a single maximum appears 
in the superconducting singlet region for Ts/U > 0.5. 

In this paper, we have clarified several characteristic 
transport properties on the assumption that the super- 
conducting gap is sufficiently large. Nevertheless we be- 
lieve that our main conclusion, i.e. the conductance max- 
imum clearly characterizes the crossover between two dis- 
tinct types of singlet ground states, holds even for gen- 
eral N-QD-S systems with a finite gap A. In this case, 
however, the crossover between two distinct types of sin- 
glet states depends also on A, as stated in §3. In near 
future, we expect that the characteristic maximum struc- 
ture found in the conductance due to the Andreev reflec- 
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tion will be observed experimentally. 
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Appendix A: Hamiltonian at |A| — > oo 

We explain how the Hamiltonian (1) can be reduced to 
an effective single-channel Hamiltonian (3) in the limit 
of |A| — > oo. For this purpose, we consider the Green 
function of the QD in the 2x2 Nambu representation, 
which takes the form 



(Tdi(t)dUt')) (Td^d^t')) 
(Td\(t)d\(t>)) (TdlmW) 



_ ( Gd],d]{t-t') Fd],d[{t - f) 

Fdl,d^(t — t') Gdl,dl(t — f) 



(Al) 



We perform the Fourier transformation of the retarded 
(advanced) Green function for eq. (A-l). From the Dyson 
equation, we have 



I. 



(A-2) 



In the interacting case, the self-energy £ r ^(e) can be 
classified into two parts 

S^( £ ) = S^ (a) ( £ ) + S r j a) ( £ ), (A3) 

where H^ a \e) is the self-energy due to the mixing be- 
tween the QD and the leads, and £^ a '(e) is due to 
the Coulomb interaction. The exact expression for Sq 
(= (Sq)^) i s written down as 



where 



-i(T N + T s l3(e)) ir S /?(e)t 
iT s (3(e)^ -i(T N + T s (3(s)) 



|e|0(| e |-|A|) , e6(\A\-\e\ 



V|A| 2 -e 2 



(A-4) 



(A-5) 



v/e 2 -|A| 2 

Here, we consider the case that the superconducting gap 
is sufficiently large, i.e. the limit of |A| — > oo. In this 
limit, taking into account /3(e) — > — 7 ===, Sg(e) is 

reduced to 



|A|^oo 



-iT N 



r s e^ s 

-iT N 



(A-6) 



We note that the off-diagonal element in eq. (A-6) can be 
regarded as a static superconducting gap induced at the 
QD, and its amplitude is given by the resonance strength 
Ts- Therefore, in the limit of |A| — ► oo the information 



about the superconducting lead can be included through 
the off-diagonal term of the Green function of the QD. 
This enables us to rewrite the Hamiltonian (1) to an ef- 
fective single-channel Hamiltonian with an extra super- 
conducting gap at the QD in eqs. (3) and (4). 

Appendix B: Bogoliubov transformation 

We perform the Bogoliubov transformation for the 
Hamiltonian (6). Using the Nambu representation, H d + 
H d c and Hn+Htn in the Hamiltonian (6) are rewritten 

as, 



N-l 

n N + h tn = ]T t„A-/ 2 + , 

(B-2) 



n=-l 



where 



(B-3) 



In eq. (B-l), £ d = e d + U/2 and f_ ltr = d a . For eqs. (B-l) 
and (B-2), we carry out the Bogoliubov transformation, 

where 



Vy,n y (-l)-!^ 



u„ 



v d u 



(B-4) 
'! ) , (B-5) 



\ u d\ 



+ N^K 1 "*)' (B ' 6) 



(B-7) 



As a result, each term of the Hamiltonian Ti^RG can 
transformed as follows, 

H° d +H^ = E d ^-yli,7_i,-l) , (B-8) 

N-l 

Un+Utn = £ ^t„A-/ 2 (ji +lalna + H.c.) 

n— — 1 a 

(B-10) 

eqs. (B-8), (B-9) and (B-10) give the Hamiltonian (10). 
We should notice that H d remains unchanged under the 
Bogoliubov transformation. 

We finally show how the Green function of the QD 
is transformed via the Bogoliubov transformation. Ap- 
plying the Bogoliubov transformation of eq. (B-4) to the 
Fourier transformation of eq. (A-l), 



G d4 {e) = U d G 1 _ ltJ _ 1 (e) U\ 



(B-ll) 



where the Green function G 1 _ in _ 1 (e) is for the Hamilto- 
nian after the Bogoliubov transformation. Therefore, we 
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end up with G 7 _ li7 _ 1 (e) without off-diagonal elements, Substituting eq. (C-5) into eq. (C-4), we have 



G 7-lt,7-lt( £ ) 

G 7-iJ.,7-iJ.( £ ) 



(B-12) 



Appendix C: Derivation of eq. (18) 

By applying the Landauer formula, we derive the linear 
conductance dl /dV\v=o- The current from the normal 
lead to the QD is given in terms of the Nambu represen- 
tation 12 - 16 

I = | J ^rv) [Gl d (e)n(e)Gl d (e) 

-2i(l-2/ JV )ImG5 id (e)] 11 (CI) 

The subscript 11 means the (1,1) element of the matrix. 
O(e) is the self-energy given by fi(e) = — S < (e) — S > (e), 
where H < ^ > \e) is the lesser (greater) self-energy for the 
QD. = f(e — eV), where /(e) is the Fermi distri- 
bution function. Similarly to S r( - a ' ) (e) in eq. (A-3), the 
self-energy Jl(e) can be classified into two parts, 



n(e) = n (e) + «£,(£), 



(C-2) 



where O (e) (Clu{e)) is due to the mixing between 
the QD and the leads (Coulomb interaction). Using 
cq. (CT), we consider the linear conductance Gy=o = 
dI/dV\v=o at zero temperature. It should be noted that 
ImS[,(0)|y = o = at zero temperature, so that the 
term including i~2[/ (e) vanishes. Namely, assuming that 
ImS[,(0)|y = o = and £lu(0)\v=o = at zero temper- 
ature, correlation effects on Gy=o at zero temperature 
enter through Rc£j} (0) | v=o ■ On the other hand, the ex- 
act formula for £lo(e) is given by 

« (e) 



-2i 



(1 - 2f N )T N + (1 - 2f s )r s ps(e) 
-(l-2/ s )r s/ 5s(e)# 



-(l-2/ s )r s/ 5s(e)# 
(1 - 2f N )T N + (1 - 2fs)T s ps(e) 



(C-3) 

where f s = f(e), f N = f(e+eV) and p s (e) = ^# 
Using the retarded (advanced) Green function in eq. 
(A-2) together with eq. (C-3), we obtain the linear con- 
ductance at the zero temperature as, 



Gv=o 



4^ 
h 



AT% \FS hdl (0)\ 2 



(C-4) 



As shown in eq. (A-l), d ^(e) is the (1,2) element of 
G d . d (e). Using C? 7 _ iTi7 _ iT (e) in eq. (14) and eqs. (B-5), 
(B-ll) and (B-12), we can rewrite d j(e) around the 
Fermi energy as 

F^ dl (s) = u d v d (G 7 _ iTi7 _ iT (e) - <5 7 _ li;7 _ a ( £ )) 



u d v d 



e-E d + iT N e + E d + iT 



N 



(C-5) 



G 



If! 4^1 
h E 2 



d (r% + Ej) 



(C-6) 



Note that the factor A^/E 2 coming from u d and v d in eq. 
(B-6) is not renormalized by the Coulomb interaction. 
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